Heavy Quark Potential at Finite Temperature Using the Holographic Correspondence 
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We revisit the calculation of a heavy quark potential in A/" = 4 supersymmetric Yang-Mills the- 
ory at finite temperature using the AdS/CFT correspondence. As is widely known, the potential 
calculated in the pioneering works of Rey et al. [^ and Brandhuber et al. 0] is zero for separation 
distances r between the quark and the anti-quark above a certain critical separation, at which the 
potential has a kink. We point out that by analytically continuing the string configurations into 
the complex plane, and using a slightly different renormalization subtraction, one obtains a smooth 
non-zero (negative definite) potential without a kink. The obtained potential also has a non-zero 
imaginary (absorptive) part for separations r > re — O.STO/ttT . Most importantly at large separa- 
tions r the real part of the potential does not exhibit the exponential Debye falloff expected from 
perturbation theory and instead falls off as a power law, proportional to 1/r^ for r > ro = 2.702/7rr. 



Heavy quark potential is a very important quantity in 
gauge theories at finite temperature. It also has a great 
phenomenological relevance in connection with experi- 
mental programs in heavy ion collisions at the Relativis- 
tic Heavy Ion Collider and at the upcoming Large Hadron 
Collider. The melting of heavy mesons in a medium is 
considered to be one of the main experimental signatures 
for Quark-Gluon Plasma formation, the ultimate goal of 
such experiments. Current analyses of available experi- 
mental data indicate that the matter formed in such colli- 
sions is strongly coupled. Thus, the study of heavy quark 
potential requires strong-coupling techniques, such as the 
Anti-de Sitter space/conformal field theory (AdS/CFT) 
correspondence [3|, U, la, |a] ■ The main goal of this work is 
to improve the current description of heavy quark poten- 
tial at finite temperature in the AdS/CFT framework. 

Until recently heavy quark potential has been cal- 
culated either analytically at small coupling using the 
perturbation theory, or numerically using lattice simu- 
lations. With the advent of AdS/CFT correspondence 
[dj y, S Ifil , it became possible to analytically calculate 
the heavy quark potential at strong coupling, albeit only 
for the A/" = 4 supersymmetric Yang-Mills (SYM) theory. 

The first calculation of a heavy quark potential in vac- 
uum for A/" = 4 SYM theory was carried out by Mal- 
dacena in [7|. Soon after [7| calculations of the heavy 
quark potential for A/" = 4 SYM theory at finite temper- 
ature appeared in [l|, |^ . In [7| the heavy quark potential 
is obtained from the expectation value of a static tempo- 
ral Wilson loop and in [l|, \^ from the correlator of two 
Polyakov loops. They are calculated by extremizing the 
world-sheet of an open string attached to the quark and 
anti-quark located at the edge of the AdS^ space in the 
background of the empty AdSs space in [7| and in the 
background of the AdSs black hole metric in [l|, ^ . 

The zero-temperature heavy quark potential obtained 
in [7| is of Coulomb type due to conformal invariance of 



TV = 4 SYM theory: 
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with A the 't Hooft coupling and cq = F^ (j) /{2n)^^^. 
Here r is the distance between the quark and the anti- 
quark in the boundary gauge theory. 

The finite-temperature heavy quark potential obtained 
in [ij, y starts out at small r being close to the vacuum 
potential of Eq. ([T]), but rises steeper than the vacuum 
potential, becoming zero at a separation r* — 0.754/7rT. 
For larger separations, i.e., for r > r* , the authors of 
[l|, 2\ argue that the string "melts" , and the dominant 
configuration corresponds to two straight strings stretch- 
ing from the quark and the anti-quark down to the black 
hole horizon. The resulting potential is thus zero for 
r > r* and has a kink (a discontinuity in its derivative) 
at r = r*. 

To clarify the definition of the heavy quark potential 
at finite temperature T let us first define the Polyakov 
loop for SU{Nc) gauge theory at spatial location f by 
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■Tr 



P exp \ ig I dr A^ (r, r) 



(2) 



with T the Euclidean time and /3 — 1/T. The connected 
correlator of two Polyakov loops can be written as [8|, ]^ 



(L(0)Lt(f))e = 






(3) 



Eq. ([3]) is the definition of singlet Vi{r) and adjoint 
Vadj (f) potentials in Euclidean time formalism. With the 
appropriate modification of Eq. ^, Eq. ^ also applies 
to AA = 4 SYM. 

To calculate the Polyakov loop correlator in AdS space 
one follows the standard prescription outlined in [l|, \^ |7| 
and connects open string(s) to the positions of Polyakov 
loops at the boundary of the AdS space in all possible 
ways. The two relevant configurations are shown in Fig.[T] 
and labeled "hanging string" and "straight strings" . 




hanging string 



straiglit strings 



FIG. 1: Two configurations of open strings corresponding to 
Polyakov loop correlator. Solid horizontal line denotes the 
boundary of the AdS space, while the dashed line denotes the 
location of the black hole horizon. 

The two string configurations shown in Fig. [T] give two 
different saddle points of the Nambu-Goto action. In the 
large- iVc large-A limit the integral over all string config- 
urations, and hence the Polyakov loop correlator as well, 
is equal to the sum of the contributions of the different 
saddle points. Actually, as was argued in [l^, the two 
straight strings on the right of Fig. [1] have Chan-Paton 
labels indicating which D3 brane each string ends on. 
Therefore the straight strings configuration actually rep- 
resents of the order of N"^ extrema corresponding to the 
different ways the two straight strings connect to Nc D3 
branes. Summing over all the saddle points we write (in 
Euclidean space) 

^hanging / O \ Qstraight 
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with S 
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NG 



and S 



straight 
NG 



the Nambu-Goto actions of 
the hanging and straight string configurations. 

Comparing Eq. (J4|) with Eq. ([3|) we conclude that the 
hanging string configuration gives Vi (r) , while the two 
straight strings stretching to the horizon give Vadjir)- 
However Vadj itself is A'^^-suppressed and repulsive, 
while Vi is of order one in A'c-counting and attrac- 
tive [8]. Renormalizing the Nambu-Goto actions in 
Eq. ^ by subtracting the actions of the string configura- 
tions at infinite quark-antiquark separations one obtains 



S 



straight 
NG.ren 



0, which implies that Vadj is zero at leading 



order in A^^. The first non-trivial contribution to Vadj is 
given by graviton exchanges between the strings in the 
bulk calculated in [10|. If exponentiated (eikonalized) , 
they would indeed give A'^-suppressed contributions in 
the exponent, as expected for Vadjir). Here we will cal- 
culate the singlet potential Vi{r). In lattice simulations 
it is usually Vi (r) which is understood as the heavy quark 
potential at finite temperature |ll| . In the real-time for- 
malism Vi (r) is given by the expectation value of a static 
(temporal) Wilson loop via 
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with the temporal extent of the Wilson loop T -^ 
c». Note that when calculating the Wilson loop (O in 
Minkowski space only the hanging string configuration 



contributes, as the quark and the anti-quark are pro- 
jected onto a color-singlet state at initial and final times. 

To find Vi{r) (henceforth referred to as V{r)) we will 
study the behavior of the hanging string solution found in 
P,[2| for r>r*. As is well-known, for r > rv = O.STO/ttT 
the string coordinates of the solution [l|, |3] become 
complex-valued. This simply indicates that the saddle 
point of the Nambu-Goto action lies in the complex string 
coordinate region: it does not invalidate the saddle point 
approximation and the results obtained with it. Simi- 
lar complex-valued solutions were recently observed by 
the authors in [iJl, where the scattering amplitude of a 
quark-anti-quark dipole on a shock wave was calculated. 
In 12] the complex- valued string coordinates were instru- 
mental for finding the unitary solution for the scattering 
cross section. Inspired by that example, below we will an- 
alytically continue the potential of [l|,l4| into the complex 
region of string coordinates. The resulting potential is 
smooth. The corresponding force on the quarks is a con- 
tinuous function of r. By modifying the ultraviolet (UV) 
subtraction we obtain a potential which is non-zero for 
all separations r. The potential develops an imaginary 
part, corresponding to the decay of the quark-anti-quark 
singlet state: similar results have been seen in finite tem- 
perature perturbation theory in 13|, |lj| • Finally, instead 
of Debye screening leading to exponential falloff of the 
potential at large distances, we find the power-law falloff 
Re[F(r)] ^ 1/r^ at large r. 

We want to calculate a temporal Wilson loop in a 
finite-temperature A/" = 4 SYM medium. We shall de- 
fine the real-time heavy quark potential in the same way 
as in [l3| . Following [l|, y] we start with the AdSs black 
hole metric in Minkowski space [3|, ll5|, llq 



ds^ = 



L' 



- l- 



dt^ + dx^ 



dz^ 



(6) 



where dx^ ~ (dx^)'^ + (dx^)^ 4" (dx^)'^, z is the coordinate 
describing the 5th dimension and L is the curvature of 
the AdSs space. The horizon of the black hole is located 
at z = Zfi with Zh = l/irT . 

We want to extremize the open string worldsheet for 
a string attached to a static quark at x^ — r/2, x"^ = 



„3 _ 



and an anti-quark at x^ — —r/2,x' 



,2 _ ^3 _ 



= 0. 



Parameterizing the static string coordinates by 

X" = [X° ^t,X^^ X, X^ = Q,X^ = 0, X^ = z{x)] (7) 
we write the Nambu-Goto action as 

SNG{r,T) = -^—T 

ZTT 

where z' — dz{x)/dx. 

The Euler-Lagrange equation corresponding to the ac- 
tion dl) is 




{2 + zz"){z^-zt)-2z'^z^ + zi)^Q 



(9) 



with z" — (Pz/dx^ . Solving Eq. ([9]) with the boundary 
conditions z{x ~ ±?'/2) = one gets 
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(10) 



where Fi is the AppeU hypergeonietric function. Here 
Zmax is the constant of integration corresponding to the 
maximum of the string coordinate along the fifth dimen- 
sion of the AdSs space, whose boundary is located at 
z — 0. It is given by the solution of the following equa- 
tion 



rco 



F 



(11) 



with F the hypergeonietric function. Indeed in the T -^ 
limit Zh ^ (X and Eq. (fTT|) gives us Zmax = rco in 
agreement with Maldacena's vacuum solution [7| . 

The action in Eq. ([5]) contains a UV divergence, which 
has to be subtracted out. Usually, the subtraction con- 
tains a finite piece as well [l|, 13, LZIj which may be 
temperature-dependent in the case at hand. Here we 
will use the following subtraction, different from the one 
used in [l|, y] : we define the quark-anti-quark potential 
by 



V{r) 



{SNG{r,T) - Re[SNGir = (X3,T)]}/r. (12) 



This subtraction insures that the real part of the poten- 
tial V{r) goes to zero at infinite separations. Our sub- 
traction p^ is consistent with that used in [7| to find 
the heavy quark potential at zero temperature. 

Using the solution from Eq. ^^ in Eqs. (O and ^2]) 
we obtain the following expression for the heavy quark 
potential of the A/" = 4 SYM theory at finite temperature: 
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Eq. P3|) below is also needed to obtain Eq. P^ . Our 
subtraction prescription resulted in the 1/ Zh term on the 
right of Eq. (fT^ instead of 2co/z;i, which would corre- 
spond to the subtraction done in [l|, I4I . Eq. ([13]) along 
with Eq. pT|) gives us the heavy quark potential as a 
function of the separation r and temperature T = I/ttz/j. 
As can be readily checked numerically, Zmax given by 
Eq. PT|) becomes complex for r > r^ = 0.870^/1, lead- 
ing to complex- valued zix) and the potential V{r). This 
led the authors of [l|, |3] to abandon their solution for 
r > re (in fact, the solution was abandoned even ear- 
lier, for r > r*). We suggest however to interpret the 
complex- valued saddle points as corresponding to quasi- 
classical configurations in the classically forbidden region 



of string coordinates. This is similar to the method of 
complex trajectories used in quasi-classical approxima- 
tions to quantum mechanics 17 1. 

The complexification of the string coordinates simply 
indicates that the saddle point of the integral over string 
coordinates becomes cornplex. According to the standard 
AdS/CFT prescription [7|, in the large- iVc large-A limit 
the integral over string coordinates is still dominated by 
the saddle point, even if it is complex. Therefore, the 
fact that string coordinates at the saddle point become 
complex does not undermine the validity of the approx- 
imation. Below we will extend the solution of Eqs. (fT3|) 
and (llip to r > Tc allowing for complex- valued Zmax- 

First we note that the extension of the solution for 
Zmax following from Eq. (|lip to r > r^ is not unique. 
Two most important roots of Eq. (|lip found in [l|, \^ are 
shown in Fig. [21 which depicts real and imaginary parts 
of Zmax as functions of separation r. There are other 
roots of Eq. ([TT|) that are not shown in Fig. [2 they are 
either negatives of the roots in Fig. [2] or the complex 
conjugates of the roots in Fig. [2] and of their negatives. 
Such extra roots lead to physically irrelevant configura- 
tions and are not shown in Fig. [2] To determine which 
root gives the correct potential from first principles one 
has to (at least) calculate quantum o(l/vA) corrections 
to the quasi-classical results shown above. While such a 
calculation is necessary, it would be rather tedious and 
is left for future work. Here we will demand that the 
correct root maps on the Maldacena vacuum solution [7| 
in the zero temperature limit. In addition we will impose 
unitarity to single out the right root. 
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FIG. 2: Real and imaginary parts of the roots of Eq. ([TT]) 
plotted as a function of quark-anti-quark separation r. 

Of the two roots shown in Fig. [2] only one (denoted by 
the solid line) maps onto Maldacena's solution behavior 
of Zmax = rco at small r [7J]. Since, on physical grounds. 



we want our potential to recover the zero temperature re- 
sult [7| at small r we will keep this root, and discard the 
other root denoted by the dash-dotted line in Fig.[21 The 
root given by the solid line develops a positive imaginary 
part for r > Vc, as shown in the bottom portion of Fig. [51 
As the complex conjugate of this root would also be a 
solution of Eq. pTj) . while mapping onto Maldacena's so- 
lution for small r, the question arises about the choice of 
one root over its complex conjugate. To select the root we 
note that the quantum-mechanical time-evolution opera- 
tor in Minkowski metric is e~*^* ~ gim[s]t_ Demanding 
that the probability of a state does not exceed one we 
obtain Im[£'] < 0, leading to Im[F(r)] < 0. This condi- 
tion allows us to single out the solid line in Fig. [2] over 
its complex conjugate as the physically relevant root. 
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FIG. 3: The real and imaginary parts of the heavy quark 
potential plotted as functions of the separation r for several 
different temperatures. We put A — 10. 

Using the solid line root of Fig. [5] in Eq. P^ . we can 
plot the real and imaginary parts of the resulting poten- 
tial. The plots are shown in Fig. [3] In the top panel of 
Fig. [3] we show the real part of the heavy quark poten- 
tial for two non-zero temperatures, along with the zero-T 
curve for comparison. One can see that the non-zero tem- 
perature curves are indeed strongly screened compared to 
the zero temperature case, but remain non-zero at all r. 
The subtraction scheme proposed above in Eq. p^ in- 
sures that Re[F(r)] given by Eq. p^ approaches zero 
at large r. The use of the subtraction scheme proposed 
in [ll y would have led to Re[V^(r)] going to a positive 
constant as r ^ oo. 

As is clear from the lower panel in Fig. [31 the heavy 



quark potential develops an imaginary part for r > Tc- 
This means the potential becomes absorptive, as the qq 
singlet state may melt in the medium. The rate of ab- 
sorption increases with r, as larger pairs are more likely 
to decay. The existence of an imaginary part in the heavy 
quark potential has been previously observed in pertur- 
bation theory in [il,[l3. While Re[t/(r)] in Fig. [3] does 
not have a kink, there is a region near r — Tc where the 
slope of the curve changes very fast. This rapid change is 
due to the potential developing an imaginary part, which 
should quickly reduce the force on the quarks. 

As one can explicitly check from Eq. (fT3)) . at small r 
we recover the zero temperature potential of [7[ 
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At large separations r we first use Eq. pT|) to write 
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Using Eq. ^^ in Eq. ^^ we obtain 



Re[U(r)] 
and 
Im[U(r)] 
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As one can see from Eq. (|T6| . instead of the exponen- 
tial falloff with r characteristic of Debye screening, which 
would have been expected from small coupling perturba- 
tion theory and which was postulated for M ~ A SYM 
theory at strong coupling in [10|, the real part of the 
heavy quark potential falls off as a power, Re[U(r)] ~ 
l/T^r^, at large r. If our hypothesis of using the com- 
plex string configurations is confirmed, this would be an 
interesting new type of screening for the potential. How- 
ever, the large negative imaginary part of the potential 
(ITT)) leads to exponential decay with time of the initial 
(built in by construction) color correlation between the 
quark and the anti-quark in the Wilson loop. 

Combining the large- and small-r asymptotics in Eqs. 
(fTi)) . P^ we can interpolate the real part of the potential 
to write an approximate formula 



Re[V{7 



V\ 



„3 
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with the new scale tq equal to 
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(19) 



Eq. (fT8|) fits the curves on the upper panel of Fig. [3] quite 
well. The parameter ro, defined by Eq. (jlSp and given in 
Eq. p^ . can be interpreted as the screening length. 

While our power-law screening Re[y(r)] ~ 1/T^r^ is 
very different from the exponential falloff due to Debye 
screening Re[F(r)] ~ exp(— tou r) in coordinate space, 
the difference is not so profound in momentum space. 
Define a Fourier transform of the potential 



V{q) = 



(fj 



■V{r) 



(20) 



Using the small-r asymptotics ()14|) one can easily show 
that at large 9 = l^l the Fourier transform of the potential 
scales as V{q) ~ l/(?^, in agreement with the standard 
perturbative result. At small q 



V{q = Q) = I (frV{r) 



1 

2^2 



^- (21) 



Such asymptotic behavior is very similar to the Debye 
screening in the infrared (IR) , given by the screened prop- 
agator Vniq) ~ 1/(q^ + "^i)) in the perturbation theory 
with Debye mass itid ^ T ^ l/j'o- Hence our V{q) 
has qualitatively the same UV and IR asymptotics as 
the standard perturbative Debye-screened potential. The 
main difference is in the shape of V{q) at finite q: unlike 
Voiq) our V{q) is concave at all q. 

To summarize, we proposed a method of calculating 
the heavy quark potential in the finite-T strongly-coupled 
A/" = 4 SYM theory in the region of separations r and/or 
temperatures T where the classical string configuration 
does not exist. We used an analogue of the complex 
trajectories method in the quasi-classical quantum me- 
chanics 17| and analytically continued the string con- 
figurations into the region of complex coordinates. This 
allowed us to obtain a potential which is physically mean- 
ingful for all values of r and T. The potential develops an 
imaginary absorptive part. We would like to stress that 
instead of Debye screening at large separations the real 
part of the potential falls off as a power of the separation, 
which is a new and never before observed phenomenon 
in relativistic quantum field theories. 
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